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Blazars show rapid and violent variabilities, which timescale are often less than a day. We stud¬ 
ied intraday variations by applying a “shot analysis” technique to Kepler monitoring of blazar 
W2R 1926-1-42 in Qnarter 14. We obtained a mean profile calculated from 195 rapid variations. 
The mean profile shows three components; one is a sharp structnre distributed within ±0.1 day of 
the peak, and two slow-varying components. This spiky-peak component reflects features of rapid 
variations directly. The profile of peak component shows an exponential rise and decay of which 
timescales are different, 0.0416 and 0.0588 day respectively. This component is too sharp to repre¬ 
sent a standard function which is often used to express blazar variations. This asymmetric profile 
at the peak is difficult to be explained by a simple variation of the Doppler factor by changing a ge¬ 
ometry of the emitting region. This result indicates that intraday variations arise from a production 
of high-energy accelerated particles in the jet. 

I. INTRODUCTION 


brightness variations in blazars are related to sizes of 
emitting regions and these speeds in relativistic jets. 
Variations, however, have a variety of timescales from 
minutes to decades. A power spectrum density (PSD) 
of blazar shows a power-law distribution, which means 
that variations of blazars show noise-like behaviors 
M- These various brightness variations in blazars 
can be happened by a variety of physical situations 
in relativistic jets. The shorter-timescale variations 
should be reflected to the physics of inner-emitting re¬ 
gions of a jet. Thus, the study of short-timescale vari¬ 
ations is important to investigate the origin of varia¬ 
tion in blazar jets. 

Blazars show rapid variations having a timescale of 
less than 1 day. These rapid variations have been re¬ 
ported in wide wavel eng ths from radio to gamma-ray 
bands; in the radio (TS|. optical Q, X-ray 0 , and 
TeV gamma-ray bands [3| . The Fermi space telescope 
has scanned the entire gamma-ray sky every hours, 
and detected a lot of large-amplitude variations as 
flares Detected flares in the gamma-ray 

band often continued for less than 1 day, and had a 
variety of shapes, not only a simple rise and decay. We 
need higher time-resolution and photon-statistics ob¬ 
servations to study the detailed feature of rapid vari¬ 
ations. 

An optical continuous monitoring of blazar 
W2R 1926±42 with a high time-sampling rate by Ke- 


Blazars have relativistic jets whose axes are closely 
directed along their line of sight an. Timescales of 


pier spacecraft [1] detected a lot of rapid variations. 
Its light curve revealed detailed shapes of numerous 
variations with large signal-to-noise ratio. We report 
general features of rapid variations having a timescale 
less than 1 day by stacking these detected variations, 
and producing a mean profile of rapid variations, so- 
called shot analysis. This paper is organized as fol¬ 
lows. Details of Kepler observation and the way of 
shot analysis are described in section 2. Observational 
results and features of the mean profile of rapid vari¬ 
ations are reported in section 3. We discuss an origin 
of rapid variations of the object in section 4, and sec¬ 
tion 5 gives several concluding remarks. 


II. OBSERVATION AND ANALYSIS 
A. Kepler data 

Kepler monitored over a hundred thousand ob¬ 
jects in Cygnus regions, and obtained continuous light 
curves with two timing settings, long (thirty-minute) 
or short (one-minute) cadences. Blazar W2R 1926±42 
is listed in Kepler target list. It has been obtained 
a continuous light curve with the long cadence since 
Quarter 11. In Quarter 14, the object had been moni¬ 
tored in the short cadence mode for 100 days. We pro¬ 
duced the calibrated “SAP_FLUX” light curve with 
one-minute time resolution by the automated Kepler 
data processing pipeline 0. 

W2R 1926±42 is classified as a low-frequency 
peaked BL Lac object at z = 0.154 estimated from 
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two absorption lines Q . Edelson et al. reported that 
there were numerous flares with timescales as short as 
a day in the Kepler light curve of Quarter 11 and 
12 [^. The PSD calculated from the light curve 
showed approximately a power-law distribution, but 
not simple. It showed a flattening at frequencies be¬ 
low 7X 10“® Hz. 


B. Shot analysis 

Frequency-domain analyses (e.g. PSD) are not easy 
to relate with physical mechanisms directly. On the 
other hand, time-domain analyses keeping phase in¬ 
formation of variations can be useful for studying 
physical mechanisms of variation. We need, however, 
large photon statistics to study variations with these 
time-domain analyses, because it is difficult to study 
detailed features of variations by using only partial 
data with observational uncertainty. Additionally, ob¬ 
served variations in blazars usually have a variety of 
shapes. Thus, it is difficult to understand general fea¬ 
tures of variations in blazars by studing only individ¬ 
ual variations. 

We apply the light curve of W2R 1926-1-42 obtained 
by Kepler to a shot analysis proposed by Negoro et al. 
to study the general features of rapid variations with¬ 
out local features of individual variations, because a 
mean profile of rapid variations calculated by the shot 
analysis can be cancelled the local features of individ¬ 
ual variations [3 ■ We analyzed following procedures 
to make a mean profile; First, we select rapid varia¬ 
tions as candidates of shots. Second, we estimate the 
observational uncertainty in the light curve. Then, we 
select the rapid variations with four times larger am¬ 
plitudes than the standard deviation of the observa¬ 
tional uncertainty after subtracting the baseline com¬ 
ponents. We define these variations as shots. 

There are two possibilities for varying the observed 
brightness, intrinsic variation of the object and vari¬ 
ation by the observational uncertainty. It is natural 
that the variation by the observational uncertainty is 
dominant rather than the intrinsic variation in shorter 
period, especially in the period between two observ¬ 
ing points which lie next to each other. We estimate 
the standard deviation a of differences between two 
neighboring points, and define a as the observational 
uncertainty, tr =17.15 count s“^. 

Rapid variations are often superposed on long-term 
variations in light curves of blazars [l^. We ap¬ 
proximate a baseline component of rapid variation by 
a second-order local polynomial fitting to the light 
curve without the period of the rapid variation. We 
subtract the calculated baseline component from the 
light curve, and extract the rapid variation. We de¬ 
tect a shot when the estimated amplitude of extracted 
rapid variation without the contribution of the base¬ 
line component is larger than our threshold, >4tT. 


Additionally, the peak time of the shot is defined at 
the time of the maximum flux among the period of 
rapid variation after subtracting the baseline compo¬ 
nent. We calculate a mean profile of detected shots 
by stacking with reference to each peak. 

III. RESULTS 

Figure [T] shows an optical light curve of the object 
obtained by Kepler. The object showed a violent vari¬ 
ability with various timescales ranging from several 
tens of minutes to over ten days, limiting for the obser¬ 
vational uncertainty. In the light curve, there are not 
only a large-amplitude long-term variation like from 
JD 2456150 to 2456160, but also a lot of flare-like vari¬ 
ations with timescale of hours. These rapid variations 
existed throughout the entire period of this monitor¬ 
ing. Figure [5] shows examples of rapid variations. We 
range these rapid variations with reference to the peak 
times of individual extrema. Figure [5] clearly shows 
that these variations had a variety of shapes. 

We detect 195 shots from the obtained light curve, 
in pursuance of the definition of shot described in sec¬ 
tion 2.2. We calculate a mean profile of these detected 
shots by appling the shot analysis. Upper panel of 
figure [3] shows the mean profile of shots without the 
data at peak time, because positive fluctuations of the 
counts at t = 0 are summed up systematically fl^ . 
There are mainly three components at the mean pro¬ 
file of shots shown in figure [31 a sharp component dis¬ 
tributed in ±0.1 day of the peak time (component 1), 
and slow-variable components ranging from —0.50 to 
—0.15 day and from 0.10 to 0.45 day (component 2 and 
3), respectively. An increase and decrease of flux in 
component 1 are approximately exponential rise and 
decay. Additionally, the profile at the peak is chang¬ 
ing from rising to decaying phases for approximately 
ten minutes. 

If shot profiles change depending on selected ampli¬ 
tudes of shots, the calculated mean profile does not 
reflect to general features of shot. We verify whether 
there is an amplitude dependence to the profile or not. 
We separate detected shots in three terms based on 
these amplitudes, 4-6 (t, 6-8cr, and over 8 ct, and cal¬ 
culate mean profiles using selected shots. Although 
calculated profiles show small differences caused by 
the limited number of shot samples, the mean pro¬ 
file of shots has no clear trend associated with these 
amplitudes. 

We estimate the systematic uncertainty of the ob¬ 
tained mean profile of shots associated with limited 
sampling by a non-parametric bootstrap approach. 
First, we resample 195 shots with replacement from 
detected shots, and calculate a mean profile from the 
resampled 195 shots. We produce 10000 pseudo mean 
profiles of shots following this procedure with differ¬ 
ent resamplings. We normalize an average of each 
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FIG. 1: Light curves obtained by iTep/er spacecraft in Quarter 14. The object monitored for 100 d with one-minute time 
resolution. 
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FIG. 2: Examples of rapid variations. These varia¬ 
tions are ranged with reference to peak times of indi- 
vidnal extrema during the period of rapid variations. 
From bottom to top, the peak times of rapid variations 
are JD 2456134.84, 2456147.23, 2456147.91, 2456151.40, 
2456152.15, and 2456153.04, respectively. 
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FIG. 3: A mean profile of detected shots. Upper panel 
shows the mean profile of shots and bottom panel shows 
the standard deviation estimated from a non-parametric 
bootstrap method. See text for detail. 


mean profile within ±1 day, and calculate standard 
deviations in each time bin. The standard deviations 
of normalized mean profiles can be regarded as the 
systematic uncertainties associated with the sampling 
of shots. The bottom panel of figure |3] shows calcu¬ 
lated standard deviations. These deviations are rang¬ 


ing from 6 to 10 count s“^. Detected components from 
1 to 3 in the mean profile of shots can be regarded as 
the real phenomena, not the artificial ones caused by 
the systematic uncertainty of the sampling of shots. 
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Component 1 reflects general features of shots di¬ 
rectly, because this component is distributed around 
the peak time. First, we apply this exponential shape 
of component 1 to a function proposed by Abdo et al.; 

F{t) = Fo + F,, (1) 


where T^. and are variation timescales of rise and 
decay phases, F^ represents a constant level underly¬ 
ing the component 1, and Fq measures the amplitude 
of the shot [l|. We evaluate its goodness of fit by 
test, = X) {Fda.ta.{ti) - Fmodei(ii)}^- The x^ of the 
best fitted function is 2278 within ±0.1 day of the 
mean profile to except contaminations of other com¬ 
ponents. On the other hand, we apply an another 
function; 


F{t) 


Fo e~OF -I- Fc, t < 0 

Fo ± Fe, t> 0, 


( 2 ) 


where Tr and Td are e-folding times of rise and decay, 
and Fc and Fq are the same in the case of function (CD- 
We also calculated the x^ of this function, x^=469. 
Figure |4] shows the applied functions with the best 
fitted parameters superposed on the mean profile and 
its residuals. Although the function ([T]) shows obvious 
residuals during the peak time in panels ”a” and “b” 
of figure m the residuals in the case of function ([2) are 
suppressed shown in panels “c” and “d”. This indi¬ 
cates that the mean profile is more spiky than the ex¬ 
pected profile from function (ED- Therefore, the good¬ 
ness of fit of the function ([2|) is more plausible than 
that of function o to represent the component 1 of 
the mean profile. 

We estimated the best-fitted parameters with a chi- 
squared test. The mean profile, however, has a sys¬ 
tematic uncertainty caused by the sampling of shots 
as mentioned above. We applied the non-parametric 
bootstrap approach to calculate the confidence level 
with the same way to estimate the errors of calcu¬ 
lated parameters. First, we calculated 10000 pseudo 
mean profiles estimated from resamples with replace¬ 
ment from detected shots. We calculated the best- 
fitted parameters against individual pseudo profiles, 
and estimated the confidence levels in these parame¬ 
ters. In table 1, we show the best-fitted parameters of 
function (I2|) to the mean profile and the ranges of 95 % 
confidence levels. We applied the Wilcoxon rank-sum 
test which was a non-parametric significance test (also 
referred to as the Mann-Whitney U-test) to the dis¬ 
tributions of Tr and Td calculated by the bootstrap 
approach [l^, [^. We confirmed the difference be¬ 
tween the Tr and Td, because the p-value was less 
than 10“^®. Therefore, the component 1 in the mean 
profile of shots has an asymmetric profile. 


TABLE I: Parameters of best-fitted function © to com¬ 
ponent 1 of the mean profile of shots 



Best value 

95% confidence level 

Tr (day) 

0.0416 

[0.0320, 0.0543] 

Td (day) 

0.0588 

[0.0399, 0.0919] 

Fo (count s“^) 

76 

[65, 88] 

Fc (count s“^) 

1508 

[1484, 1537] 
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FIG. 4: Best-fitted functions snperposed on mean profiles 
of shots. Panels “a” and “c” show a mean profile of shots 
and the best-fitted functions 0 and 0. Panels “b” and 
“d” show residnals between the mean profiles and the es¬ 
timated best-fitted fnnctions. 


IV. DISCUSSION 


We obtained the optical continuous light curve of 
blazar W2R 19426±42 with one-minute time resolu¬ 
tion by Kepler spacecraft. The object showed violent 
variability and a lot of rapid variations with timescales 
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less than a day. We detected 195 rapid variations 
as shots of which amplitude were larger than 4cr af¬ 
ter subtracting these baseline components, and ap¬ 
plied to the shot analysis. The mean profile produced 
from detected shots shows three components, one fast- 
spiky (component 1) and two slow-varying compo¬ 
nents. Component 1 shows an asymmetric profile, 
faster-rise and slower-decay features with the spiky 
but smooth-connected peak. 

It is poorly understood whether rapid variations 
are intrinsic phenomena or apparent one caused by 
a geometrical changing in the jet. There are several 
models that flux variations are explained as apparent 
brightness variations, for example varyin g th e Doppler 
factor for changing the viewing angle [24| or grav¬ 
itational lensing effect These models, however, 
expect that the averaged variation profile is almost 
symmetric in a simple situation, because the Doppler 
factor should be changed symmetrically in the aver¬ 
aged variation. In other words, the rise and decay 
timescales should be equal. Estimated rise and decay 
timescales of rapid variations, however, are different 
described in section 3. Thus, these models can be 
ruled out in the case of rapid variations. Thus, rapid 
variations may come from authentic phenomena. It is 
plausible that there is a particle acceleration during 
the rapid variation, and then, higher-energy particles 
increase in the emitting region of the variation. 

The synchrotron cooling timescale Tgyn is repre¬ 
sented as, sec, where 

B is a strength of magnetic field, E is an observed 
energy band, and 5 is the Doppler factor [2l|, |2^. If 
the dissipation of high-energy particles in rapid vari¬ 
ations is caused by the synchrotron cooling, Tsyn in 
the rest frame can be represented using the decaying 
timescale T^, Tgyn = 5 Td/{l-\-z). We estimate 5 of 5.8 


from observed Td of the mean profile of shots, where 
E is 2.25 eV and assuming B as 0.5 G which is typical 
value among the gamma-ray detected BL Lac objects 
m- The mean profile reflects common features of 
rapid variations. Thus, the estimated S should be a 
typical value of inner regions where rapid variations 
happen. 


V. CONCLUSION 

The optical continuous light curve with one-minute 
time sampling obtained by Kepler revealed that the 
mean profile of rapid variations almost showed ex¬ 
ponential rise and decay. Rise and decay timescales 
of shot profile, however, are different, and the profile 
shows asymmetric profile. A particle acceleration pro¬ 
cess can produce this asymmetric variations. There 
are several scenarios which can explain the particle- 
acceleration mechanism causing the rapid variations; 
shock- in-j et scenario [l^ . magnetic reconnection sce¬ 
nario [Q. The shot analysis is also feasible to 
study the spectral feature of variations, because of 
large signal-to-noise ratio. Unfortunately, Kepler per¬ 
formed only one-band monitoring. Spectral and fur¬ 
ther observational studies are needed to completely 
understand the mechanism of rapid variations. 
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